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Abstract — In this work, a MIMO detection problem is con- 
sidered. At first, we derive the Generalized Likelihood Ratio 
Test (GLRT) for arbitrary transmitted signals and arbitrary 
time-correlation of the disturbance. Then, we investigate design 
criteria for the transmitted waveforms in both power-unlimited 
and power-limited systems and we study the interplay among 
the rank of the optimized code matrix, the number of transmit 
diversity paths and the amount of energy integrated along each 
path. The results show that increasing the rank of the code matrix 
allows generating a larger number of diversity paths at the price 
of reducing the average signal-to-clutter level along each path. 

I. Introduction 

We consider a detection system wherein both the receiver 
and the transmitter, which are not necessarily co-located, are 
equipped with multiple widely spaced antennas. The scenario is 
that typical of a Multiple-Input Multiple-Output (MIMO) radar 
architecture, as outlined in [1] and analyzed in [2], [3], but may 
also describe the tasks of a sensor network with widely spaced 
nodes, wherein both the transmit and the receive array elements 
exhibit no coherence. In this situation, the MIMO architecture 
can be advocated to improve the detectability of targets with 
fluctuating radar cross section either through an increase of the 
diversity order [1], [2] or through beam-forming in the signal 
space [3] or through a combination of these two strategies: the 
degree of freedom which allows jumping from one strategy 
to another is the structure of the transmitted signals, and in 
particular the rank of the space-time code matrix used at the 
transmitter 

In the above context, this paper makes the following contri- 
butions. At the receiver-design level, we first derive the Gener- 
alized Likelihood Ratio Test (GLRT) for arbitrary transmitted 
waveforms and arbitrary time-correlation of the disturbance. 
Building upon [2], [3], two design criteria for the space-time 
code are then presented and discussed. The figure of merit is the 
Mutual Information (MI) [4] between the received signals and 
the channel vectors generated by a prospective target, whose 
maximization is related to the ability to estimate unknown 
characteristics of the target as shown in [5], and also to the 
maximization of the detection probability {Pd) for a fixed false 
alarm probability (P/a) as pointed out in [2], [3]. We examine 
both the case of power-unlimited systems, wherein a constraint 
is forced upon the received average signal-to-disturbance ratio, 
and the case of power-limited systems, wherein the constraint is 
instead on the transmitted energy. Finally, closed-form formulas 
for Pfa and Pd are given, showing that there is an inherent 
trade-off between number of diversity paths and amount of 
energy integration along each path, which is a direct confirma- 
tion of the well-known fact that no uniformly optimum (e.g., 
for any signal-to-disturbance ratio) coding strategy exists. The 
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Fig. 1. Considered scenario. 



tool that can be used for trading diversity order for integration 
is the code-matrix structure, and particularly its rank. 

This paper is organized as follows. In Section|ll] the problem 
under investigation is described. In Section|lIIl the GLRT-based 
detector is derived. Finally, the diversity-integration trade-offs 
for different code design criteria are analyzed in Section HV] 

II. Problem Statement 

Consider a detection architecture where the transmitter and 
the receiver consist of M and L widely spaced antennas, 
respectively, and are not necessarily co-located as depicted in 
Figure[T] We assume the following general model for the signal 
transmitted by the m-th antenna (node) 

N 

Sm{t) = an,mK{t) , < i < , m = 1, . . . , M , (1) 

where N is the signal-space dimension, {4>n{i)}n=i is an 
orthonormal 

(ai,m, ■ ■ ■ ,aN 



basis of the considered signal space, am 
^ is the codeword associated with the m- 
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transmitter, ||am|| is the energy transmitted by antenna m and 
Ts is the time duration of the transmitted signal. Subsumed by 
the above model is the case of a pulsed-MIMO radar, wherein 
4>n{i) = P{t -{n- l)T), with p{-) a pulse of duration Tp <T 
and T = Ts/N the pulse-repetition time (PRT) [2]. 

The signal observed at the L receive nodes may contain or 
not the echo of a target scattering located at a given distance. 
The bandwidth of the transmitted signal induces a partition of 
the controlled area in a finite number of range cells, directly tied 
to the delay with which a target echo is heard at the receiving 
nodes. Denoting by rg{t) the signal received at the ^— th receive 
node at a given delay tq, we thus have: 



Hi 



M N 

EE 

m— 1 n— 1 



{t - To)a,n,e + Wi{t) 



(2) 



Ho : wt{t) 



for a. = 1, . . . ,L. In the above equation, which subsumes 
the model outlined in [2], we{t) is a Gaussian process which 
models the overall disturbance originating from the super- 
position of receiver noise and spurious scattering from the 
surrounding environment (denoted as clutter) which typically 
exhibits time correlation. In writing down dU, the following 
three assumptions have been made, a) The transmitted signals 
are narrowband [1], namely all of the receive sensors see the 
target as belonging to the same range celfl b) The spacing of 
the antennas at both the transmitter and the receiver is wide 
enough as to allow angle diversity [1], [6], i.e., to generate 
different aspect angles which justifies the use of the LM 
different coupling coefficients a„i.e for m = l,...,Af and 
£ — 1, . . . , L. c) The target is either stationary or has a known 
Doppler shift which is compensated for at the receiver [7]. 

In principle, Analog-to-Digital (A/D) conversion would re- 
quire undertaking a Kelly-Root expansion of the received 
(vector) signal [8], which would inevitably lead to an A/D 
stage depending on the clutter covarianc^ This would in turn 
lead to a hardly implementable receiver and would prevent 
any subsequent step towards adaptive systems, wherein such 
a covariance is assumed unknown [3]. As a consequence, 
we adopt the customary - albeit sub-optimum - approach of 
projecting the received observations onto the signal space, 
which yields the following binary hypothesis testing problem: 
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(3) 



where ctg = {au, . . . , a^.i)^ is the Af— dimensional vector 
modeling the scattering of the target - hit by the M transmitted 
waveforms in ([U - towards the th receive node for £ = 
1, . . . , L; A ~ (ai, . . . , a^/ ) is the N x M code matrix 
containing the A^— dimensional codeword transmitted by the 
TO— th transmit node as its m— th column; tr{AA^} is the 
total transmitted energy; finally, wi represents the overall 
disturbance at the th receive antenna and is modeled as 
an A^— dimensional Gaussian vector with known (full-rank) 
correlation matrix M. 

For future reference, we define A = min{A^, M} and 5 ~ 
rank{A}, with 1 < 5 < A. Also, we define the received 
average signal-to-clutter ratio (SCR) as 



SCR = 



1 
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:^tr{A^M-iA}, 
(4) 

where the expectation is over the realizations of the vector a = 
vec{(ai, . . . , chl)}, tr{-} denotes trace and the assumption has 
been made that E[afQ;|^] = (j'^Im for £ = 1, . . . , L. 

Given the above ingredients, we first derive a distribution- 
free test for (|3]l, namely a test whose implementation does not 
require prior knowledge of the target parameters: in this work, 
we assume that detection takes place at a fusion center wherein 

'This assumption can be relaxed througii cell synchronization. 

^This becomes the conventional Karhunen-Loewe expansion to be under- 
taken at each receive antenna when, as assumed throughout this study, the 
clutter exhibits temporal coiTelation, but no spatial correlation. 
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Fig. 2. GLRT statistic computation. 

un-quantized versions of the r^'s are made available. Next, 
we investigate two design criteria for the code matrix A, so 
as to elicit the interplay between the transmitted waveforms 
si(i), . . . , SM{t) and detection performance. 

III. Detector Structure 

Lacking prior information as to the target statistics, the 
Neyman-Pearson test cannot be implemented, whereby a wise 
design strategy for the detector is the Generalized Likelihood 
Ratio Test (GLRT) [7] which, under the model Q, yields: 



max In 

ai,...,aL 



/(ri, . . .,rL\Hi,ai, . . ..oll) 



/(ri,...,ri|77o) 
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where /(ri, . . . , |i?i, ai, . . . , ccl) and /(ri, . . . , rL|i?o) 
denote the conditional densities of the observations under the 
two alternative hypotheses, while ?/ is the detection threshold, 
to be set based upon the desired false alarm probability, Pfa 
say. The maximization involved by (|5]l is obviously tied to 
the structure of A, which can be in principle arbitrary, and 
in particular to its rank i5, which in turn may take on any 
integer value from 1 to A. Straightforward derivations lead 
to the following general form of the GLRT, which holds for 
arbitrary code-matrix A: 
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wherein Pm-^/^a represents the orthogonal projector onto the 
range span of the matrix M^^^^A. 

Let BDC^ be the singular value decomposition (SVD) [9] 
of M^^/^A, where B is a unitary N x N matrix containing 
the left singular vectors, D is an N x M matrix containing 
the singular values on the main diagonal with [DJi i > • • • > 
[-D]a,a > 0, and C is a unitary M x M matrix containing 
the right singular vectors. Under rank-(5 coding, the orthogonal 
projector in (|6]l can be written as Pm-^/^a — BsBf , wherein 
Bs is the N x 5 matrix containing the first 5 columns of B. 
Three relevant cases are subsumed by this solution. 

1) Under rank-1 coding, we have that A = aq^ with 
a A^— dimensional and q unit-energy Af— dimensional; 
hence, in this case the orthogonal projector reduces to 



M 



2) Under rank-A (i.e., fuU-rank) coding and < M, we 
have Pm-i/2a = In- 



3) Finally, under rank-A coding and N > M, we have 

Notice that M^^/'^r^ represents the output of a filter aimed 
at whitening the disturbance impinging on the ^— th transmit 
antenna, whereby the computation of the GLRT test statistic in 
(|6]l may be undertaken based on the scheme of Figure|2] 

The performances of the above test can be expressed in 
closed - albeit implicit - form under very general conditions. 
In particular, under rank-(5 coding, we have: 
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where P — ol^ A Acxi and Qm{-, •) denotes the 

generalized Marcum Q function of order m. Notice that ^ 
and ^ generalize previous results derived in [2] for 5 = A 
and in [3] for S — 1, respectively. 

IV. Code Design and Performance Bounds 

Code optimization requires determining a meaningful per- 
formance measure to be optimized. Paralleling the arguments 
of [2], we assign a prior to the vectors of the scattering 
coefficients; in particular, we assume that a consists of LM 
independent and identically distributed circularly-symmetric 
Gaussian random variables with equal variance cr^. At this 
point, a reasonable design criterion in the GLRT frame- 
work is to choose A so as to maximize the MI between 
the received observations ri, . . . , tl and the channel vectors 
ai, . . . ,aL under Hi, i.e., /(ri, . . . , tl; ai, . . . ,aL\Hi) = 
iylog2 [det(JAr + alM^iAA^M-^)] [4]. Indeed, since 
the GLRT amounts to writing the conditional likelihood given 
the unknown target parameters and replacing them through 
their maximum-likelihood estimates performed based upon 
the available observations, maximizing the MI between the 
received signals and the channel vectors is expected to improve 
the GLRT performance [5]. 

To make the maximization non-trivial, a set of physical 
constraints on A should be added. In the following we discuss 
two relevant cases: I) a power-unlimited system, wherein we 
constrain the normalized received average SCR in (IDl; II) a 
power-limited system, wherein we constrain the transmitted 
energy per signal dimension, i.e., ti{AA^}/N. In both cases, 
an additional constraint can be imposed upon the rank of A. 



A. Power-unlimited system 

The problem to be solved is the following: 

det (Jat + alM-iAA"M-^ 



max log2 

A 



s.t. SCR^ ^Ti-{A"M-^A} < 
rank{A} < 



(9) 



with 1/ > and 6* e {1,2,..., A}. Let BAB" be the 
spectral decomposition of the positive semi-definite matrix T = 



M-^AA^M^K where A = DD" = diag{Ai, . . . , Ajv} 
with Ai > ... > Aa? > 0. Since the rank constraint forces 
Ag+i = . . . = Xn = 0, the problem (|9]l can also be recast as 



2 

S.t. SCR = ^ y A, < 1/ 



(10) 



A straightforward application of the Jensen inequality [4] 
allows now deriving the following condition for optimality: 



A, 



iVSCR 



J = l, 



(11) 



which in turn implies that the optimal A should be such that 



(12) 



According to ( fTTl i and (fT2] l. the optimal space-time code must 
have rank{A} = 6 = 9 (i.e., the code must exploit the 
maximum possible number of degrees of freedom as permitted 
by the rank constraint) and SCR = ly; whereby, both inequality 
constraints in (|9]l can equivalently be replaced with a strict 
equality. The following remarks are now in order. 

For a target SCR and a given rank constraint 9, the above 
coding strategy amounts to generating 9 independent and iden- 
tically distributed diversity paths at each receive antenna; also, 
each path enjoys the same average signal-to-clutter ratio given 
by NSCR/9. Notice that increasing the number of diversity 
paths 9 at each receive antenna (i.e., increasing the rank of A) 
comes at the price of reducing the received average SCR per 
path. This intuition will be farther exploited in what follows. 

The solution in (fTTT i coincides with the solution that we 
would obtain if, under a definite rank constraint, we adopted the 
lower Chernoff-bound on the detection probability as objective 
function. This result was first observed in [2] for full-rank 
coding (i.e., 9 = A). 

Finally, we point out that the optimal code matrix complying 
with ([T2I is not unique. In fact, let Ua^aV^ be the SVD 
of A, where Ua is a unitary N x N matrix containing the 
left singular vectors, is an TV x M matrix containing the 
singular values on the main diagonal and Va is a unitary 
M X M matrix containing the right singular vectors. Also, let 
Um^'m^m 1'^ '■h^ spectral decomposition of M^^, where 
Um = {um^i, ■ ■ ■ ,um,n) is an N X N unitary matrix and 
All = diag{l/AA/,i, . . . , 1/Aa/,jv} with 1/Aa/,i > ... > 
^/Xm,n > 0. Then, (fT2l l can be rewritten as 



UMA-y^U^UA-^A 



'NSCR 
~9ar 



Bo 



(13) 



It is clear that the unitary matrix Va can be arbitrarily chosen 
since it does not come into play in (T3[ . Also, it is instructive 
to consider the class of optimal solutions obtained by setting 
Ua — Um- In this case, the right singular vectors of A 
are matched to the N eigenvectors of the noise matrix M 



which in turn define as many orthogonal modes in the signal 
space. The position of the 9 non-zero singular values in 
determines which subset of orthogonal modes is employed for 
transmission. In particular, let n(l), . . . ,7i(6') be the row (or 
column) indexes corresponding to the non-zero singular values 
in S^. Condition (fT3] l is fulfilled by setting 



Power-unlimited: = 4, M 



■'A\n{])Mj) 
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(14) 



It is seen from ( fT4b that more energy must be allocated to more 
interfered modes in order to equalize their received average 
SCR's. Among all the possible choices, setting n{i) = j for 
j — 1, . . . ,9 minimizes the total transmit energy tr{S^S^} 
required to achieve a given SCR since transmission takes place 
along the least interfered modes. 

1) Analysis and trade-offs: The detection performance of 
the GLRT under rank-^ optimal coding in (fTTT i has a simple 
closed-form expression. Indeed, since the constraint on SCR 
allows generating exactly 9L independent diversity paths with 
one and the same received average SCR given by A^SCR/6', 
the test statistic (|6]l is a Gamma random variable with shape 
parameter 9L and scale parameter 1 + A^SCR/6' under Hi, 
implying that 



Prl 



exp 



l + iVSCR/6' 



eL-i , 

y If ^^'i 

^ l\ Vl + A^SCR/^y 



The inherent trade-off between number of generated diversity 
paths BL and amount of energy integration granted on each 
path A^SCR/6' is visible already at this stage, and is confirmed 
by the asymptotic behavior (for increasingly large SCR) of the 
detection probability, i.e.: 
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(15) 



Notice that no strategy is uniformly superior in terms of 
detection performance, but the optimal value of 9 depends 
upon the operating SCR, as also confirmed by Figure[3] Indeed, 
while for large SCR the asymptotic behavior in ( fTsT i indicates 
that maximizing the diversity order (i.e., choosing 6* = A) 
amounts to maximizing the detection probability, in the low 
signal-to-clutter regime this trend does not hold anymore. To 
have an intuitive justification of this fact, recall that increasing 
the number of diversity paths eventually leads to more and 
more constrained target amplitude fluctuations. On the other 
hand, ([8]l reveals that the detection probability is in the form 
Pd = [P(/3,?7)] = [P(/3rms/3o,??)], where /3™^ is the 
root mean square value of the random variable /3, while /3o 
is a random variable with unit rms value. ry), regarded 
as a function of /3, is a cumulative distribution function and, 
therefore, exhibits a sigmoidal shape (i.e., is U— convex in the 
region /Jrms 0, while being n-convex in the region firms 
oo). Under these circumstances, more and more constrained 
fluctuations (i.e., lower dispersion ratios of the parameter f3) 
typically result in larger values of P^j in the n— convexity 
region, while being detrimental in the U— convexity region. 
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Fig. 3. Probability of miss versus SCR for the optimal coding in lilt . 

B. Power-limited system 

In power-limited systems, the problem to be solved is 



max 

A 



l0g2 



dct / 



s.t. tT{AA"} < N8 
rank{A} < 9 



(16) 



where £ > is the available transmit energy per signal 
dimension and 9 E {1,2,..., A}. A solution to (fTSI l was 
derived in [3] for 9 ~ N < M. Also, maximizing the objective 
function in ( fTSI l for 9 — 1 amounts to maximizing SCR in 
(lUi and the solution to this problem under a transmit energy 
constraint can again be found in [3]. In the following, the 
solution to ( fTSI l for arbitrary values of 9 is discussed. 
Notice first that the problem ( fT6b can be recast as 



max 

T>-0 



log2 [det {In + <tIT)] 



(17) 



s.t. ti{M^TM2} < N£ 
rank{T} < 9 

where we recall T = BAB^ . Since the rank constraint forces 
Xe+i = . . . = Xn — 0, ([TtI i is also equivalent to 



max 
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l0g2(l+^aA,) 



(18) 



S.t. tr{BeAeBf M} < N£ 



where Ag = diag{Ai, . . . , Ae}, Bq — (6i, . . . , bg) and bj is 
the j— th column of the unitary matrix B = {hi,..., bM)- 
Observe that the choice of Bg only affects the feasible region 
of the problem ( fTSl l. As a consequence, the optimal Bg must 
ensure the largest feasible set for Ag which rules the value 
of the objective function. This implies that the first 9 columns 
of B must be equal to the eigenvectors um.i, ■ ■ ■ ,UM.e of 
M corresponding to the 9 smallest eigenvalues [10, Lemma 
9.H.l.h]. The last N — 9 columns of B can be arbitrarily 
chosen (provided that B is unitary) since they do not come 
into play in the optimization ( fTSl l. Therefore, we can assume 
in the following that B = Um, which in turn implies 



Ua = Um and T = L/mAj,/ P^A^/ C/iJ with Pa = 
Tia'^a — diag{pA4, ■ • ■ ,PA,e, . . . , 0}. Hence, (fTsT l reduces 
to the following standard water-filling problem [4] 
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PA.i>0,...,PA.f)>0^ \ 



s.t. 



PA.j < N£ 



(19) 



whose solution is simply given by 



PA.j = max I 0, — 



Y,PA,j^N8. (20) 




Some remarks are now in order. Under tr{ AA^} < N£ pig. 4. Probability of miss versus 7 for the optimal coding in (20), 



and rankjA} < 9, the optimal code matrix A — Ua'^aVa 
is recovered (up to a right multiplication by an AI x A/ 
unitary matrix) by forcing Ua = U]\i, [EA]j.j — y/PA.j for 
j = 1,...,9 and pA.j given by ( |20l i. and [S^]j,j = for 
j = + 1, . . . , A. This solution (which subsumes the results 
in [3] for the special cases 9 = 1 and 9 = N < M) has 
a nice physical interpretation. The optimized code discards 
the noisiest N — 9 modes defined by M in the signal space. 
Also, since the remaining 9 orthogonal directions present a 
different disturbance level, more energy is allocated to more 
reliable modes according to ( |20] i. The number of activated 
modes depends not only upon the upper rank constraint 9, but 
also upon the available transmit energy £, the target strength 
(7^ and the eigenvalues \m,i < ■ ■ • < ^M,e of the disturbance 
covariance matrix. If £a'^/\M,e 00, then paa > ••• > 
PA,e > and S = rank{A} = 9. Instead, if not enough energy 
is available for transmission or if the target is weak, additional 
modes may be switched off. Also, notice that the activated 
modes have in general different received average signal-to- 
clutter ratio's given by a-'^pA.j / Xm.j for j = 1,...,S < 9 
(more on this infra). 

We emphasize that the situation here is dramatically different 
from what observed in the previous section under a received av- 
erage SCR constraint, wherein the transmitter always activates 
the largest possible number of orthogonal modes permitted by 
the rank constraint and pumps into them as much energy as to 
make these paths equivalent at the receiver side. 

1) Analysis and trade-offs: In Figure|4] we report 1 — 
(probability of miss) for the optimal water-filling solution in 
(EQIi. We assume Pfa = lO^^, M = N = A and L = 1; also, 
we consider an exponentially-shaped disturbance covariance 
matrix with [AfJ^j = o-^O.?'*^^', whose eigenvalues are 
{0.203, 0.318, 0.754, 2.725}aj. To elicit the joint effects of the 
transmit energy constraint and the target strength, the curves 
are plotted versus the parameter 7 = ^cr^/cr^. 

Notice first that, if 7 ^ 0, the optimal code matrix acti- 
vates only the least interfered mode, independent of the rank 
constraint 9; in this case, the test statistic (|6j is a Gamma 
random variable with shape parameter L and scale parameter 



1 + (j'^NE / Xm,i under Hi, implying that 



Pd = exp <^ - 



1 + alNS/XM 



NS/X 



M,l 



This result confirms the basic intuition that, when the transmit 
energy is scarce or if the target is weak, we should give up 
diversity and concetrate all the energy along the least interfered 
direction in the signal space in order to maximize the received 
average signal-to-disturbance ratio. A similar result also holds 
in MIMO point-to-point wireless channels [6]. 

On the other hand, if 7 is large enough, relaxing the rank 
constraint 9 allows generating S > 1 orthogonal modefl In 
our example, if —6.1 dB < 7 < 2.3 dB up to three modes can 
be activated; instead, full rank coding becomes possible only if 
7 > 2.3 dB. For a fixed £, it is clear that increasing the number 
6 of modes (i.e., diversity branches) again cames at the price 
of reducing the energy transmitted (and therefore the average 
SCR ratio received) along each path. As a consequence, there 
is no strategy which is uniformly optimal for all values of 7. 
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